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Abstract. Solar supergranulation plays an important role in generating and 
structuring the solar magnetic field and as a mechanism responsible for the 11- 
year solar cycle. It is clearly detected within SOHO/MDI Dopplergrams, from 
which a variety of properties may be derived. Techniques that extract spatial, 
temporal and kinematic characteristics and provide comparisons for the two 
most recent solar minima are described. Although supergranule lifetimes are 
comparable between these minima, their sizes maybe slightly smaller during the 
recent minimum. 



1. Introduction 

Supergranules are the largest scale component of solar convection readily seen 
with present observing techniques. They are typically ^30 Mm across a nd have 
been observed to live from anywhere between 1-2 days ( Duvalll Il980l ). Their 



relationship with local magnetic fields have been seen via diffusion studies and 
spatial observations of both magneto grams and CallK ch romospheric data. 

Since its inception, SOHO/MDI ( Scherrer et aLlfTgOSl ) Dynamics Runs have 



provided yearly 60-day sets of 1-minute Dopplergrams. Analysis of the 1996 
data has quantified a wide variety of supe rgranule characteristics (for example, 
iHathawav eralll2000l : iBeck fc Sch^l200(]l 1. 



The most recent solar minimum has sparked much interest due to its lengthy 
minimum. Many studies are currently underway to understand the causes and 
consequences of this recent minimum. 

This paper extends the analysis of supergranule characteristics to 2008 data 
and the results are compared to similar studies performed on 1996 data. 



2. Methods 

A series of data reduction processes were applied to the Dopplergrams to isolate 
the Doppler velocity signals of supergranules. The analysis was then split into 
studies of supergranule velocities, cell sizes and their 1/e lifetimes. 



2.1. Data Reduction 

The p-mode oscillations were removed by performing a weighted average over 
31 1-minute de-rotated Dopplergrams ( Hathaway 19881 ). The resulting time- 
series was sampled every 15 mi nutes, producing 96 im ages per day. Next, the 
global flow fields were removed (|Hathawavlll987l . Il992l ^. The remaining Doppler 
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signals are due to the photospheric convective velocity field. These convective 
Dopplergrams were remapped to heliographic coordinates followed by a pro- 
.j ection onto the spher ical harmonics to produce their respective power spectra 



( Hathaway et al. 2000l ). Instrumental calibration artifacts were also removed. 



2.2. Supergranule Velocities 

Hathawav et al.l ( 20021 ) performed a study of supergranule velocity flows using 



the 1996 data and compared the relative strength of the radial and horizontal 
components. That study was extended to investigate both the 1996 and 2008 
data, detailed by Williams & Pesnell (2010, in preparation). 



2.3. Supergranule Lifetimes 

Correlation techniques have been used to follow supergranule patterns in the 
Dopplergrams as they traverse across the solar disk and measure the change in 
the pattern over the period of a day. A data strip, extracted near the equator, 
was selected from an image. A similar strip window in a subsequent image, at a 
given time-lag from the original, was moved longitudinally prograde across the 
image and the correlation noted for each pixel shift. This process was performed 
for all images within a single day, resulting in a two-dimensional array (pixel- 
shift vs. time-lag) of correlation coefficients. The region of decaying positive 
correlation pertaining to the evolving supergranule pattern was isolated using a 
Radon transform, from which the maximum correlation coefficient, C(r), at a 
given time lag, r, was found. A plot of In C(t) vs. r was linearly fitted, the slope 
of which provided the decay coefficient of the correlation (Figure [1]) . The inverse 
of this decay coefficient gave the 1/e decay-time of the supergranule pattern. 

Performing this process for the 1996 & 2008 data resulted in a series of 
decay-times for each day. Treating each decay-time as an independent result, a 
harmonic mean provided the average decay-time for each data-set. 



2.4. Supergranule Sizes 

Using the method described bv lHathawav etah ( 200d ) . power spectra were pro- 
duced for each Dopplergram contained in the 1996 and 2008 datasets. Mean 
spectra were produced for both years by averaging over all the individual spec- 
tra. A best-fit to the supergranule feature {I < i < 250) was deduced, using 
a modified Lorentzian function (Figure [2D . Although the fit deviates from the 
data at higher wavenumbers, this is of no consequence to the study as only 
the behavior near the peak is of interest. The location of the peak in the fit- 
ting determines Imax and the FWHM of the supergranule distribution. From 
the peak and FWHM, the typical diameter and size range of supergranules was 
determined. 



3. Results 



The supergranule velocity analysis was performed out to an angular distance of 
p ~ 30 degrees. The 2008 data produced radial and horizontal RMS velocity 
components of Vr = 40 m/s and Vh = 303 m/s, with Vr/vh = 0.13. The 1996 
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Figure 1. A log-linear plot (crosses) of correlation coefficient vs. time-lag 
represents the change of a supergranule pattern during a single day. A linear 

fit to the plot (solid line) extracts the decay coefficient for the convection 
pattern and, subscqticntly, its 1/e decay-time. The squares signify the points 
between which the fit was performed. For this particular day (in this case 
May 4, 2008), the decay coefficient is 0.0643, giving a decay-time of 16 hours. 
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Figure 2. A 60-day average of MDI convection spectra from 2008 (points) 
is fit with a modified Lorentzian. The fit (dashed line) is produced within a 
given range (1 < ^ < 250). The function is extended throughout the full data 
range for illustration purposes (dotted line). The FWHM (solid line) is shown 
within the supergranule feature. 



data produced Vr = 24 and Vh = 278 m/s, with Vr/vh = 0.09. These values have 
an error of =b 1 m/s. 

For 2008, the correlation analysis gave a mean decay coefficient of 0.056 it 
0.001 correlation coefficient units per hour. Calculating decay-times from each 
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coefficient, tlie harmonic mean calculation gave a 1/e decay time of 18 liours. 
Tliis value was also derived using the same process for the 1996 data. 

For the supergranule size analysis, the 2008 dataset gave a peak wavenumber 
of £ = 124 ± 1 and an FWHM of At = 168 ± 5. These results respectively 
correspond to an average supergranule diameter of 35.2 Mm and a diameter 
range of 18.7-67.2 Mm. For 1996, the peak wavenumber was found to be ^ = 
121 lb 1 and an FWHM of Ai = 160 it 3, corresponding to an average diameter 
of 36.1 Mm with a diameter range of 19.5-68.3 Mm. 



4. Analysis, Discussion and Conclusion 

While the 1/e lifetimes are found to be the same for both years, the supergranule 
cell sizes tend to be slightly smaller during 2008. A significant discrepancy, 

however, is found within the velocity analysis. 

Although the 1996 results are in line with those found bv lHathawav et"al 



( 20021 ) ■ th e 2008 velocity values a re considerably greater. Image defocusing dur 



ing 1996 ( Korzennik et all |2004| ) may influence the velocity values contained 



within the Dopplergrams by altering the resolution of the instrument. Active 
regions present on the disk rnay co ntribute to spurious velocity values within 
the calculation ( Liu &: Norton|[200ll ). Any images displaying significant activ- 



ity are currently removed from the analysis, whereas future analyses will use 
magnetograms to mask out these regions so the images may be included. 

Comparing the past two solar minima, the decay rate and decay-times of 
the supergranule pattern have not changed. Missing images within a day's data 
currently means removing that day from the analysis, which reduces the sam- 
pling for the averaging. Future work will update the analysis process so that 
loss of data within a day will not mean that the whole day needs to be skipped. 

The size analysis is much less problematic although it can be improved using 
a statistical analysis of the extracted parameters for each Dopplergram-derived 
spectrum. This is planned for the future along with producing similar statistics 
for the lifetimes and velocity analyses. 
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